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1. Introduction 

In previous work [[TJ, 0, some of us outlined the concept for a magnetic field monitoring system 
for a neutron electric dipole moment experiment in which the spatial dependence of the vector 
components of the magnetic field (and, hence, the gradients dBj/dxj) within the inaccessible in¬ 
terior measurement region are reconstructed solely from exterior measurements. In this paper, we 
show results from first prototyping tests of such a system. The results highlight the potential for 
the implementation of an improved system in an upcoming neutron electric dipole moment experi¬ 
ment to be carried out at the Spallation Neutron Source at Oak Ridge National Laboratory (i.e., the 
SNS nEDM experiment [0 0]). The work presented here complements techniques that have been 
developed to extract interior gradients from exterior measurements in scalar magnetometers 0 . 


2. Methodology 


Detailed discussions of the concepts underlying our field monitoring system were published pre¬ 
viously El-0- In brief, the basic idea is as follows. In a region of space containing no sources 
of current or magnetization (key assumptions underlying our concept), the magnetic scalar poten¬ 
tial obeys the Laplace equation V 2 4>v/(x) = 0, for which the solution can be written in spherical 
coordinates according to the well-known multipole expansion: 


+(. 


0,0) = ^ rP™ (cos d ) ai m cos (in(f))+be m sin (m(j>) 

£—0 in —0 


(2.1) 


Here, the P"'(cos6 ) denote the associated Legendre polynomials, and a/ m and b( m denote a pri¬ 
ori unknown expansion coefficients. The magnetic field can then, of course, be written as B = 
—V4>m(x). Thus, if values for the magnetic field components B t can be determined at a number of 
exterior points, such measurements can be employed to determine values for the expansion coeffi¬ 
cients ai m and b( m , which in turn permits reconstruction of the scalar potential in the interior region 
(up to an arbitrary constant). Taking the gradient of the expansion equation ( |2.1| ) then uniquely 
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Table 1 . Sensitivity of the field components (B x ,B y ,B z ) to the ap m and bp m expansion coefficients in terms 
of their associated basis functions, f( m j{x) and h( ln j(x), in the expansion of the magnetic scalar potential. 
Some trivial basis functions are not listed, e.g., fooj, whose associated basis function in the scalar potential 
expression is a constant and thus does not appear in the expressions for the B,. Also, expansion coefficients 
associated with basis functions which are linearly dependent on lower-order terms are redundant and do 
not need to be determined. Note that a ‘0’ indicates that particular component offers no sensitivity to that 
particular expansion coefficient. 
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determines the spatial dependence of the vector field components everywhere within the interior 
region. 

Translating from spherical to rectangular coordinates, we write the components Bi{x) in terms 
of the ci£ m and b( m expansion coefficients as 

Bi (*) = Yj a( 'n fom, i (*) + b im *hn, i {*), i<E{x,y,z}. (2.2) 

where we use fe m .j{x) and ly m j(x) to denote the basis functions in our expansion. Explicit expres¬ 
sions for f(:m.i(x) and hp m j(x) written in terms of ( x,y,z ) rectangular coordinates (i.e., the coordinate 
system most compatible with typical experimental field measurements) are given in Appendix B of 
Ref. [ 0 . 

Table [I] summarizes the sensitivity of (B x . B y ,B z ) to the ap m and be m expansion coefficients (up 
to £ = 3 ) in terms of these />,„,■ and h( mi basis functions. The point here is that measurements 
of (B x .By.B z ) provide (mostly) redundant information on the ap m and bp m expansion coefficients, 
as most of these coefficients are common to all three (B x ,B y ,B z ) components; however, as can be 
seen in the table, certain coefficients can only be determined from measurements of a particular 
component (e.g., b\ \ can only be determined from a measurement of By). 

If one desires reconstruction of all the (B x ,B y ,B z ) field components to a certain order (writ¬ 
ten in terms of the ap m and bp m coefficients up to that order), it is possible via measurements of 
any of the (B x ,B y ,B z ) components providing appropriate sensitivity to these coefficients, with the 
number of required measurements equal to the number of ap m and bp m coefficients up to that order. 
As a specific example, suppose one desired to reconstruct (B x ,B y ,B z ) to order (£.m) = ( 3 , 3 ) in 
the scalar potential. As can be inferred from the table, such a reconstruction would require the 
determination of only 12 coefficients, ci \o through <233, and, hence, only 12 exterior measurements. 
Five of these coefficients (020, «30, «3i, ^31, (>32) could be determined from a measurement of 
any of (B x ,B y ,B z ) at any location in space where their respective basis functions are non-zero, four 
of these (<221, ^21, ^22, ^33) could be determined from any two of (B x ,B y ,B z ), and three of these 
(a 10, an, (hi) would require a measurement of a specific component (e.g., a 10 is sensitive only to 
B z ). Thus, one can optimize the choice of external measurements according to the needs of the par- 
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Figure 1 . Left panel: A top view of the prototype array design. Right panel: Schematic figure of the 
experimental apparatus. A scaled version of the fiducial volume is also shown. 


ticular experiment, with a careful choice of such providing for maximal information in the (Lm) 
reconstruction from a minimal number of exterior measurements. 

However, if the experiment is somehow constrained such that exterior measurements of only 
one of the field components is possible, then one can only reconstruct that particular component. 
For example, as can be seen in Table [TJ, B x does not provide sensitivity to aio, bn, and l> 2 \ , coeffi¬ 
cients which would be needed for the reconstruction of B y and B-. Continuing this example, with 
12 exterior measurements of B x , one could reconstruct B x up to a somewhat higher (£,m) = (4,1) 
order, but would sacrifice the ability to reconstruct B y and B z . 


3. Experimental Apparatus and Procedure 

A prototype field monitor system consisting of twelve single-axis fluxgate magnetometer probes 
was deployed within the magnetic field environment of an optimized cos 6 coil surrounded by 
multiple layers of magnetic shielding [0,0| developed as a part of prototyping studies for the SNS 
nEDM experiment. 

The single-axis fluxgate magnetometer probes, obtained from Stefan-Mayer Instruments 0, 
were mounted on a cylindrical-like support structure consisting of four aluminum rods, as shown in 
figure HJ, with the cylindrical axis of this support structure oriented along the axis of the cos 0 coil. 
For the purposes of this first demonstration, the arrangement of the probes was rather simplistic; 
all of the probes were oriented along the cos 9 coil field direction (i.e., perpendicular to the axis of 
the cos 9 coil), with four probes mounted (identically) on each of three planes oriented, as shown 
in figure |T], perpendicular to the coil axis. 
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To set the geometric scale for these first tests, note that the radius of this cos 0 coil is 32.4 cm, 
with a length of 214.6 cm [J7p. The effective radius of the support structure for the field monitor 
probes was 12.85 cm, with the three planes separated by a distance of 30.5 cm. During data 
taking, the cos 6 coil was energized to a nominal field of ~ 220 nrGauss at the coil center. The 
readings in the twelve field monitor probes were then recorded simultaneously at one point in time. 
These twelve readings of the field component along the primary field direction, which we define 
to be B x , were then fitted to the multipole expansion in equation ( |2.2| ) to determine the at,,, and 
bf m coefficients up to order (4,1), as discussed in the example given in the previous section. The 
magnetic field in the region interior to the prototype array was then calculated using these fitted 
coefficients. The accuracy of the reconstruction was then assessed by comparing the results of 
the reconstruction with explicit measurements of the field in the interior region earned out using a 
triple-axis fluxgate magnetometer probe m mounted to an automated magnetic mapper system. 
Note that the readings in the twelve field monitor probes were obtained prior to the measurements 
with the automated mapper system; we did not later average the monitor probes’ readings over 
some period of time. 

As is well known, a fluxgate magnetometer probe may read a non-zero value in a true “zero 
field” environment due to an offset associated with its geometry and its electronics. We extracted 
the offsets of each of our single-axis fluxgate magnetometer probes (relative to that of the triple-axis 
probe on the automated mapper system) via a procedure in which a large number (of order ~ 700) 
of readings from the mapper probe (conducted with the cos 6 coil de-energized) obtained in the 
interior region of the field monitor probes’ support structure were fitted to a high-order multipole 
expansion. These fits were then extrapolated to the single-axis probes’ locations, and any such 
differences between the extrapolated fitted values and the single-axis probes’ actual readings were 
attributed to an offset. 

4. Results 

Results from our reconstruction of the interior magnetic field component B x (i.e., the component 
oriented along the cos 6 coil’s primary field direction) along the three axes (7-ax is along the axis of 
the cos 0 coil’s cylindrical support structure) are compared with direct measurements of the field 
obtained with the field mapper system in figure Note that the reconstructed and measured field 
components B x were normalized to a value of 1.0 at the coil center. As can be seen there, the 
agreement is good, with the fractional gradient along the x-axis, (dB x /dx)/B x , on the level of 10 4 
cm -1 . Note that whereas the reconstructed fields were based on the readings obtained in the 12 
single-axis probes at some instant in time, the measured fields (with the single triple-axis probe) 
were obtained over a time scale of order ~ 1 hour (due to the necessity of having to move the 
triple-axis probe to the different locations). Thus, the comparison between the reconstructed and 
measured fields may be degraded in the presence of time-varying background fields (which are 
not completely shielded), either on short time scales (e.g., between point-to-point measurements of 
the interior fields by the triple-axis probe), or on long time scales over the duration of the interior 
field measurements. To minimize the possibility of any such issues, the data were taken at night 
(data taking during the day was nearly impossible as a result of daytime operations of an overhead 
crane and delivery trucks) and the triple-axis probe was moved point-to-point by the automated 
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Figure 2. Comparison between the field component B x (i.e., along primary field direction) reconstructed 
from the exterior measurements in the field monitor array (dashed lines) and direct measurements obtained 
with the automated mapper system (dots). The three panels show results along the x-, y-, and z-axes. (Note 
that the error bars on the data points are smaller than the symbol size.) 


mapper system in as expeditious manner as possible. Nevertheless, as can be seen in figure 0, some 
anomalies were observed, such as that at y = — 3 cm. 

Averaged over a volume corresponding to the SNS nEDM experiment’s measurement volume 
(i.e., over a dense grid of points uniformly tilling this volume, of which a subset are the data 
points along the central axes shown in figure @), the reconstructed and measured fields agreed 
to better than ~ 1.0%. The reconstructed fractional gradients, ( dB x /dx)/B x , ( dB x /dy)/B x , and 
(■ dB x /dz ) /B x , averaged over this half-scale version of the measurement volume were —1.1 x 10 4 
cm' 1 , —1.8 x 10^ 4 cm' 1 , and —0.3 x 10~ 4 cm' 1 , respectively, which are to be compared to the 
direct measurements, which were —1.4 x 10~ 4 cm" 1 , —1.2 x 10^ 4 cm' 1 , and —0.1 x 10 -4 cm' 1 , 
again showing good agreement. 

Improved results for the SNS nEDM experiment will ultimately be required. There, the goal 
will be to monitor the field gradients to the level of 10 5 cm' 1 or better. Although this first proto¬ 
type demonstrating our method does not yet meet this criterion, we are currently investigating the 
feasibility of employing vector field probes with smaller noise and offset characteristics, and in a 
next-version prototype we will also carry out a full optimization of the probe locations for maximal 
sensitivity to the successive higher-order (l,m) terms in the nrultipole expansion. 
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